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ABSTRACT 


The  orientation  of  a  spacecraft  when  it  is  in  a  stable  equilibrium 
state  is  studied.  The  spacecraft  contains  a  single  momentum  wheel  which 
stores  angular  momentum.  The  axis  of  the  wheel  is  not  necessarily  parallel 
to  a  principal  axis  of  the  spacecraft.  For  an  arbitrary  but  fixed  speed  of 
the  wheel  relative  to  the  spacecraft  body,  it  is  shown  that  there  may  be  one, 
two,  or  three  stable  equilibrium  points.  At  each  of  these  points,  the  body 
may  spin  about  an  axis  which  is  fixed  in  both  spacecraft  and  inertial  coor¬ 
dinates.  The  orientation  of  this  axis  in  spacecraft  coordinates  can  be 
determined  from  expressions  in  this  report.  A  threshold  can  be  deter¬ 
mined  such  that,  if  the  speed  of  the  wheel  relative  to  the  spacecraft  is 
larger  than  this  threshold,  there  is  only  one  stable  equilibrium  point.  This 
information  can  be  used  to  determine  the  eventual  orientation  of  the  space¬ 
craft  provided  it  is  designed  such  that,  if  it  becomes  seriously  misoriented, 
the  momentum  wheel  drive  motor  automatically  holds  the  wheel  speed  to  a 
predetermined  value. 
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STABLE  EQUILIBRIUM  ORIENTATION  OF  A  SPACECRAFT 


WHICH  CONTAINS  A  MOMENTUM  WHEEL 


1  .  lntr  od  uction 

The  orientation  of  a  spacecraft  when  it  is  in  a  stable  equilibrium  state 
is  studied  in  this  report.  The  spacecraft  is  assumed  to  contain  a  single 
momentum  wheel  with  its  spin  axis  in  an  arbitrary  direction  and  a  damper  for 
dissipation  of  energy.  Two  situations  are  already  well  understood. 

(a)  The  wheel  is  spinning  at  such  a  speed  that  the  magnitude  of 
its  angular  momentum  equals  that  of  the  entire  spacecraft.  In  its  unique  stable 
equilibrium  state,  the  spacecraft  will  orient  itself  such  that  the  momentum 
vector  of  the  wheel  coincides  with  the  momentum  vector  of  the  entire  space¬ 
craft  which  is  fixed  in  inertial  coordinates. 

(b)  The  wheel  has  stopped  spinning  relative  to  the  spacecraft  body. 
There  are  two  stable  equilibrium  states.  In  each,  the  axis  of  maximum  moment 
of  inertia  of  the  spacecraft  is  oriented  in  one  of  two  possible  directions  along 
the  momentum  vector  and  the  spacecraft  is  spinning  about  this  axis. 

In  each  of  these  two  situations,  when  the  spacecraft  is  not  in  a  stable 
equilibrium  state,  the  damper  dissipates  energy  and  drives  the  spacecraft 
toward  a  stable  equilibrium  state.  *  In  this  stable  equilibrium  state,  the  damper 
dissipates  no  energy  since  the  spacecraft  is  either  stationary  or  spinning  about 
a  fixed  axis. 

Intermediate  situations  between  these  two  extreme  situations  are 
analyzed  in  this  report.  The  momentum  wheel  speed  relative  to  the  body  is 
assumed  to  be  held  constant  by  a  drive  motor  and  the  orientation  of  the  space¬ 
craft  when  it  is  in  a  stable  equilibrium  state  is  investigated.  It  will  be  shown 
that  for  an  arbitrary  but  fixed  wheel  speed  in  any  stable  equilibrium  state  the 

*  In  theory,  no  dissipation  takes  place  when  the  spacecraft  is  in  an  unstable 
equilibrium  state,  but  this  case  is  ignored  in  this  report  because  in  practice 
a  disturbance  torque  will  perturb  the  spacecraft  from  an  unstable  equilibrium 
state. 
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body  will  rotate  about  some  axis  which  orients  itself  along  the  spacecraft 
momentum  vector.  The  damper  dissipates  no  energy  while  the  spacecraft 
spins  about  this  axis.  This  axis  is  not  necessarily  the  wheel  axis  or  a  princi¬ 
pal  axis  of  the  spacecraft.  If  the  wheel  speed  is  larger  than  some  threshold, 
it  will  be  shown  that  there  is  only  one  stable  equilibrium  state  and  the  orienta¬ 
tion  of  the  spacecraft  when  in  this  state  can  be  determined.  If  the  spacecraft 
is  designed  so  that  when  it  becomes  misoriented  the  momentum  wheel  is  auto¬ 
matically  driven  at  a  predetermined  speed,  the  methods  of  this  report  can  be 
used  to  determine  the  equilibrium  orientation  of  the  spacecraft. 

The  precise  evaluation  of  the  expressions  of  this  report  requires  a 
knowledge  of  h,  the  magnitude  of  the  spacecraft  angular  momentum.  For 
example,  the  threshold  varies  linearly  with  h.  However,  in  many  cases,  the 
direction  of  the  axis  in  spacecraft  coordinates  about  which  the  body  spins  is 
insensitive  to  fairly  large  variations  in  h  as  long  as  the  wheel  speed  remains 
above  the  threshold.  Therefore,  the  results  of  this  report  are  useful  when  h 
is  known  only  approximately. 

When  two  of  the  principal  moments  of  inertia  are  equal,  the  results  of 
this  report  require  some  special  interpretation  since  some  of  the  expressions 
become  infinite.  The  best  way  to  handle  this  situation  is  to  let  the  moments  of 
inertia  differ  by  a  small  quantity  and  observe  what  happens  when  this  quantity 
approaches  zero.  In  practice,  no  two  moments  of  inertia  are  exactly  equal, 
so  the  results  for  equal  moments  of  inertia  are  not  emphasized  in  this  report. 

For  the  special  case  where  the  spacecraft  body  is  symmetric  about  the 
wheel  axis,  the  threshold  value  of  the  wheel  speed  is  determined  in  Ref.  Z 
using  a  somewhat  different  approach. 

Z.  Formulation  of  the  Problem 

Consider  a  spacecraft  which  is  rigid  except  that  it  contains  a  momentum 
wheel  that  is  free  to  rotate  about  a  shaft  parallel  to  an  axis  Z  in  spacecraft 
coordinates  which  passes  through  the  spacecraft  center  of  mass.  Define  the 
direction  of  positive  Z  such  that  a  vector  representing  angular  velocity  of  the 
wheel  with  respect  to  the  rest  of  the  spacecraft  is  in  the  direction  of  positive  Z. 
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The  momentum  wheel  is  assumed  to  have  symmetry  and  maximum  moment  of 
inertia  about  its  spin  axis,  and  all  parts  of  the  spacecraft  other  than  the  wheel 
will  be  referred  to  as  the  body.  Two  additional  coordinates  X  and  Y  are 
chosen  so  that  XYZ  forms  a  right-hand  orthogonal  coordinate  system  fixed 
to  the  spacecraft  body  and  centered  at  the  spacecraft  center  of  mass.  If  the 
momentum  wheel  were  clamped  in  any  position,  the  rigid  spacecraft  would  have 
three  principal  axes.  Denote  these  axes  x,  y,  and  z,  which  form  a  right-hand 
orthogonal  coordinate  system  such  that  the  positive  Z  axis  is  in  the  octant 
x  0,  y  ^0,  and  z  £  0,  and  such  that  I  ^1  and  I  ^  I  where  1,1,  and  I 

X  y  x  z  X  y  z 

are  the  principal  moments  of  inertia  of  the  spacecraft  (if  the  wheel  is  clamped) 
See  Fig.  1  . 

Let  the  moment  of  inertia  of  the  wheel  about  its  axes  be  C  and  the 

angular  velocity  of  the  wheel  with  respect  to  the  body  be  a  vector  oo  along  Z. 

The  spacecraft  has  a  total  angular  momentum  h  in  some  arbitrary  direction. 

The  vector  h  can  be  broken  into  two  components,  a  vector  Coo  along  Z  and  a 

vector  h  — Coo.  Let  the  components  of  Coo  along  x,  y,  and  z  have  magnitude 

a  Coo,  a  Coo,  and  a  Coo  where  oo  =  oo  and  let  the  components  of  h  have  magni- 
x  y  z  -  z  ?  l  ~  l  £ 

tude  b  h,  b  h,  and  b  h  where  h  =  I  h  I.  Then  a  +  a  +  a  =l,b  +b  + 

^  x  y  z  —  x  y  z  x  Y 

b  =  1 ,  and,  because  of  the  way  x,  y,  and  z  were  defined,  a  £  0,  a  >  0, 
z  77  x  y 

and  a  ^  0. 
z 

3.  Points  of  Minimum  Kinetic  Energy 

In  this  section,  we  consider  a  problem  where  a,  a,  a,  I,  I,  I, 

x  y  z  x  y  z 

C,  h,  and  oo  are  fixed  and  examine  the  minima  of  the  function 


2T(b  ,  b  ,  b  )/hZ  =  (b  -oa  )Z/l  +(b  -oa  )Z/l  +  (b  -0a  )Z/l  +;Z/C 
x  y  z  '  x  x  x  y  y  y  z  z  z 

u) 


where  T  is  the  kinetic  energy  of  the  spacecraft  and  o  =  (-jp).  Equation  1  is 
a  standard  representation  of  the  kinetic  energy  in  terms  of  components  of 
angular  momentum  and  principal  moments  of  inertia.  See  Ref.  1  or  3. 


For  each  point  on  the  unit  sphere,  b^  +  b^"  +  b^  =  1,  the  value  of  T 
can  be  computed  using  Eq.  1.  The  assumption  that  oo  is  fixed  implies  that  a 
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motor  drives  the  wheel  at  a  constant  speed  relative  to  the  rest  of  the  space¬ 
craft.  As  suming  that  energy  is  dissipated  in  the  body  by  the  damper,  regardless 

of  the  initial  values  of  b  ,  b  ,  and  b  (initial  orientation  of  h  in  the  spacecraft 

x  y  z  —  r 

coordinate  system),  the  spacecraft  will  eventually  orient  itself  so  that  b  ,  b  ,  and 

x  y 

b  take  on  the  values  of  the  coordinates  of  a  minimum  of  T.  In  this  stable  equili¬ 
brium  state,  the  direction  of  h  is  fixed  in  both  inertial  coordinates  and  in  the 
spacecraft  coordinates.  Therefore,  the  spacecraft  body  is  either  stationary 
with  respect  to  an  inertial  frame  of  reference  or  rotating  about  h.  In  this 
report  the  terms  minimum  and  stable  equilibrium  point  are  used  interchangea¬ 
bly.  At  a  minimum,  the  damper  dissipates  no  energy  since  the  body  is  either 
stationary  or  spinning  about  a  fixed  axis. 

There  are  several  situations  in  which  knowledge  of  the  components  of 
h,  or  equivalently  the  direction  of  the  axis  about  which  the  body  spins,  in 
spacecraft  coordinates  is  valuable.  Two  situations  are: 

(a)  The  spacecraft  loses  power  for  a  while  (e.  g.  ,  while  in  the 
earth1  s  shadow)  so  that  the  momentum  wheel  stops  spinning  relative  to  the 
body.  Power  is  regained  and  the  momentum  wheel  is  automatically  accelera¬ 
ted  to  a  predetermined  speed  to  relative  to  the  body.  In  many  practical  cases, 
the  magnitude  of  the  angular  momentum,  h,  will  be  much  larger  than  the 
magnitude  of  angular  momentum  changes  caused  by  disturbance  torques  during 
the  power  outage  so  that  the  direction  and  magnitude  of  h  remain  approxi¬ 
mately  fixed  in  inertial  space.  Knowledge  of  the  components  of  h  in  the  space¬ 
craft  coordinate  system  would  be  extremely  helpful  in  bringing  about  complete 
recovery  of  the  spacecraft  from  misorientation  suffered  while  the  momentum 
wheel  was  without  power. 

(b)  Due  to  faulty  ejection  of  the  spacecraft  from  the  launch 
vehicle  or  large  disturbance  torques  caused  by  failure  of  gas  jets  in  the  open 
position,  the  magnitude  and  orientation  of  h  in  inertial  coordinates  is  unknown 
but  the  wheel  is  automatically  held  by  the  motor  to  a  predetermined  speed  go 
relative  to  the  body.  Knowledge  of  the  orientation  of  h  in  spacecraft  coordi¬ 
nates  would  be  helpful  as  a  first  step  in  determining  the  orientation  of  h  in 
inertial  space  and  eventual  spacecraft  recovery. 
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It  is  shown  in  the  appendix  that  the  function  T  (b^,  b  ,  b^)  always  has 

one  and  only  one  minimum  in  the  region  b^  >  0,  b^  >  0,  b^  ^  0.  T  may  have 

zero,  one,  or  two  minima  in  the  region  b  <0,  b  ^  0,  b  >0  and  there  are 

x  y  z 

no  minima  in  any  other  regions.  Two  sufficient  conditions  that  T  have  a 
unique  minimum  are 


2  .coC.2  ^ 

0  =  (-T-)  > 


(I  -I  ) 

x  y 


”  '  '  Ha  ..I. >2/5  +  (a.  .I..)2/3‘  3 


x  y 


y  x 


j 


and 


(^)2  >  — 


(I  -I  ) 

X  z 


h  2/3 


L' 


(a  I)  +  (a 


X  z 


i  3 

1  X  J 


(2) 


(3) 


If  either  Eq.  Z  or  3  is  satisfied,  T  will  have  one  and  only  one  minimum.  Let 
t  =  C/h  times  the  positive  square  root  of  the  smaller  of  the  two  expressions  on 
the  right-hand  side  of  Eqs.  Z  and  3.  The  coordinates  of  this  minimum  satisfy  the 
relation  givenin  Table  1  for  various  ranges  of  the  parameters  a^,  a  ,  and  a^. 

Table  1 


Conditions  on  a  ,  a  ,  a 
x  y  z 


Coordinates  of  Unique  Minimum  Satisfy 


a  = 

i, 

a 

=  a 

=  0* 

b  = 

i , 

b 

=  b 

=  0 

X 

y 

z 

X 

y 

z 

a  = 

i. 

a 

=  a 

=  0* 

b  = 

l, 

b 

=  b 

=  0 

y 

X 

z 

y 

X 

z 

a  = 

i, 

a 

=  a 

=  0* 

b  = 

l, 

b 

=  b 

=  0 

z 

X 

y 

z 

X 

y 

a  = 

o, 

a 

>  o, 

a  >  0* 

b  = 

o, 

b 

=  Da 

I  b  /tea  I  +(I  -I  )b  1 

X 

y 

z 

X 

y 

y 

z  z  z  y  z  y  z 

b  = 

+ 

l  - 

b  2 

y 

\ 

z 

a  =0,  a  >0,  a  >0 
y  x  z 


b  =  0 ,  b  =  oa  I  b  /fpa  I  +(I  -I  )b  1  , 
y  z  zxx  xz  xzx 

b  =  +  /l  -b  2 

X  *\  z 


*  For  these  four  cases,  a  necessary  condition  for  T  to  have  a  unique  minimum 
is  that  either  Eq.  2  or  3  with  >  replaced  by  >  be  satisfied. 
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Table  1  (Cont.  ) 


a  =0, 

a  >  0,  a  >0 

b 

=  o, 

b  =oa  I  b  /  [  pa  I  +  ( I 

z 

x  y 

z 

y  y  x  x  x  y  : 

b 

=  +  J 

1  -b  2 

X 

y 

a  >0, 

a  >  0,  a  >0 

b 

=  oa 

I  b  /[oa  I  +(I  -I  )b  1 

X 

y  z 

z 

z 

XX  X  Z  X  z  x 1 

b 

=  oa 

I  b  /[oa  I  +(I  -I  )b  1 

y 

y 

xx  x  y  x  y  x 

b  =  +  i _b  2-b  2 

X  %  y  z 


Given  the  values  of  the  spacecraft  parameters,  a  threshold  t  of  wheel 
speed  can  be  determined  from  Eqs.  2  and  3  such  that  if  to  >  t  the  vector  h 
has  only  one  stable  equilibrium  point  in  the  spacecraft  coordinate  system.  In 
this  equilibrium  position,  the  components  of  h  can  be  determined  from  the 
relations  of  Table  1.  If  go  <  t,  there  may  be  as  many  as  three  equilibrium 
positions.  In  each  equilibrium  position,  the  spacecraft  body  may  rotate  about 
h.  These  positions  may  be  studied  using  the  expressions  of  the  appendix.  It 
is  desirable  that  the  spacecraft  be  designed  to  keep  u)  >  t  so  that  the  additional 
equilibrium  positions  will  not  exist. 

4.  Example:  Wheel  Axis  Along  a  Principal  Axis 

Consider  a  spacecraft  composed  of  two  rigid  bodies.  See  Fig.  2.  The 
first  body  contains  a  momentum  wheel  with  its  spin  axis  parallel  to  Z.  The 
second  body  rotates  very  slowly  (e.  g.  ,  once  per  day  for  spacecraft  in  synchro¬ 
nous  equatorial  orbit)  about  Z  with  respect  to  the  first  and  Z  is  a  principal 
axis  of  both  bodies.  This  slow  rotation  is  caused  by  a  drive  which  for  long 
periods  of  time  (e.  g.  ,  twenty-four  hours)  does  not  require  any  sensor  input 
from  outside  the  spacecraft  so  that  during  any  particular  short  time  interval 
of  interest  the  spacecraft  can  be  considered  to  be  a  rigid  body  with  one  prin¬ 
cipal  axis  along  axis  Z.  Assume  that  the  moments  of  inertia  normalized  by 
dividing  by  the  constant  I^  are 

Iv/l„  =  5.  330  +  0.  038  cos  26, 

A  Z 

I  /i  =  5.  160  -  0.  038  cos  20, 

Y  Zj 
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a  nd 


0.  38  sin  Z0 


■XY7' 


z 


where  0  is  the  angle  of  body  1  relative  to  body  Z. 

Consider  the  satellite  position  0=0.  In  this  case,  the  X,  Y,  and  Z 
axes  are  the  principal  axes, 


a  nd 


I  /i  =  5.  368, 
X  z 


I  /i  =5.  1ZZ. 

y  z 


I  he  case  described  here  is  a  =  1,  a  =  a  =  0.  Equations  Z  and  3  imply  that, 
u)C  /  z  x  y 

il  0  =  (-T— )  >  1— (I  /I  )  =  0.814,  a  unique  minimum  occurs  at  b  =  1,  b  = 
h  z  x  z  y 

b  =  0.  Therefore,  if  the  speed  of  the  wheel  is  always  kept  high  enough  so 
that  8  1.4$  of  the  angular  momentum  of  the  spacecraft  is  due  to  the  wheel 
spinning  relative  to  the  body,  only  one  stable  equilibrium  position  exists.  In 
this  position,  the  body  may  be  spinning  around  the  Z  axis. 


If  c  <  0.  814,  the  results  of  the  appendix  show  that  there  may  be  two 
equilibrium  positions  at 


b  =  0, 

y 

b  =  ol  /(I  -I  )  =  d/0.  814, 

z  X  X  z 

and  _ 

b  =  ±  1  -(o/0.  814)Z  . 

X  \ 

Consider  the  satellite  position  0  =  tt/Z.  Equations  Z  and  3  imply  that, 

p 

if  o  =  (-T— )  >  0.  811,  a  unique  minimum  occurs  at  b  =  1,  b  =  b  =0.  At  this 
h  z  y  x 

stable  equilibrium  position  the  body  may  spin  about  the  Z  axis. 

5.  Example:  Wheel  Axis  Which  may  not  be  Along  a  Principal  Axis 

Consider  a  spacecraft  composed  of  two  rigid  bodies.  See  Fig.  3.  The 
first  body  contains  a  momentum  wheel  with  its  spin  axis  parallel  to  Z.  The 
second  body  rotates  with  respect  to  the  first  very  slowly  about  Y  and  Y  is  a 
principal  axis  of  both  bodies.  With  respect  to  this  rotation,  the  situation  is 
the  same  as  that  of  the  previous  example  so  that  for  stability  analysis  the 
spacecraft  can  be  considered  to  be  a  single  rigid  body  which  contains  a 
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momentum  wheel.  Assume  that  the  moments  of  inertia,  normalized  by 
dividing  by  the  constant  Iy,  are 

Ix/ly  =  7.  4  +  0.  089  cos  Z6, 

Iz/ly  =  6.  89  -  0.  089  cos  29, 

and 

1^.^/ly  =  0.089  sin  20 


where  9  is  the  angle  between  the  two  bodies  which  is  considered  to  be  a  con¬ 
stant  at  any  particular  time.  The  coordinates  x,  y,  and  z  are  chosen  so  that 
z  is  along  Y  and  x  is  in  the  XZ  plane  displaced  an  angle  cp  from  X.  The 

condition  I  =0  determines  c f>.  By  substituting  coordinate  transformations 

xy  7 

into  the  definition  for  the  moments  of  inertia,  we  obtain  the  relations 


and 


V1  z  +  !x*  +  lz  !xz 


I 

X 


sin  cp 


^x  lz  !xz 


0, 


which  can  be  solved  for  I  /i  ,  I  /i  ,  and  cp. 

y  z  x  z 

Consider  the  case  0  =  n/ 4.  Then 


I  /i  =6.  875 
y  z 

I  /i  =7.  415 
x  z 

and 

cp  =  9.6°. 


Therefore,  a  =  0,  a  =0.  167,  a  =  0.  985.  Equations  2  and  3  imply  that,  if 
C  z  x  y 

p  =  ( )  >  0.  05,  there  is  only  one  point  of  minimum  energy.  At  this  point, 


b  =0  and  b  and  b  satisfy 
x  y  x  7 


=  6.  35  obx/(D  +  0.  47  b^) 


(4) 


and 


b*  =  <5» 
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This  follows  from  Table  1.  Therefore,  if  the  wheel  speed  is  kept  high  enough 
so  that  more  than  5#  of  the  angular  momentum  of  the  spacecraft  is  due  to  the 
wheel  spinning  relative  to  the  body,  the  momentum  vector  has  only  one  stable 
equilibrium  point  in  spacecraft  coordinates.  In  this  equilibrium  position,  the 
vector  h  is  not  parallel  to  the  wheel  axis.  The  body  may  spin  about  an  axis 
along  h. 


If  0  <  0.  051,  there  may  be  one,  two,  or  three  equilibrium  points 

depending  upon  whether  Eqs.  4  and  5  have  one,  two,  or  three  solutions  in 

the  region  b  =0,  l^b  ^0,  1  >  b  >  - 1 .  This  follows  from  the  appendix, 

z  y  x 

Consider  the  case  0=0.  Then 


I  /i  =  7.  489, 

X  z 

I  /i  =6.  801, 

y  z 

and  cp  =  0°  so  that  a  =  1  and  a  =  a  =0.  If  n  >  0.  092,  there  is  only  one 

y  x  z 

point  of  minimum  energy  located  at  b  =1,  b^  =  b  =  0.  In  this  stable  equi¬ 
librium  state,  the  body  may  spin  about  the  Z  axis. 


If  p  <  0.  09Z,  there  are  two  stable  equilibrium  points  at  b^  =  0,  b^  = 
o/0.  092,  and  =  i^l-fc/O.  092)2. 

Consider  the  case  9  =  n  /2.  Th  en 


I  /i  =  7.  31  1 , 

X  z 

I  /i  =  6.  979, 
y  z 

and  CD  =  0°  so  that  a  =  1  and  a  =  a  =0. 

y  x  z 

point  of  minimum  energy  located  at  b^  =  1, 
librium  state,  the  body  may  spin  about  the 
two  stable  equilibrium  points  at  b^  =  0,  b^ 


If  o  >  0.  0454,  there  is  only  one 

b  =  b  =0.  In  this  stable  equi- 
x  z  1 

Z  axis.  If  p  <  0.  0454,  there  are 
=  o/O.  0454  and  b  =  ±Jl  -  (o/o.0454)2. 


6.  Conclusions 


For  t  <  u)  <  h/C  where  t  is  some  threshold,  u)  is  the  momentum  wheel 
speed,  C  is  the  wheel  inertia,  and  h  is  the  magnitude  of  the  spacecraft  angu¬ 
lar  momentum,  there  is  only  one  equilibrium  state  of  the  spacecraft.  If  the 
momentum  wheel  axis  is  not  parallel  to  a  principal  axis  of  the  spacecraft,  in 
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the  equilibrium  state  the  spacecraft  body  spins  about  an  axis  which  is  not 
parallel  to  either  the  wheel  axis  or  a  principal  axis  when  t  <  u>  <  h/C.  If  the 
wheel  axis  is  parallel  to  a  principal  axis,  in  the  equilibrium  state  the  body 
spins  about  an  axis  which  is  parallel  to  the  wheel  axis  when  t  <  u>  <  h/C.  The 
value  of  t  and  the  orientation  of  the  spin  axis  of  the  spacecraft  body  in  the 
equilibrium  state  can  be  determined  using  the  expressions  derived  in  this 
report. 
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Fig.  2 
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Fig.  3 
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APPENDIX 


o  a  -  0, 
x 


1.  All  minima  of  T(b^,  b^,  b^)  satisfy  the  following  expressions.  If 


b  =  0 
x 


oa  I 

^  x  y 


(A  -  ]  > 


oa  I 

by  =  (I“-IX) 

x  y 


-  1  <  b  ^  1 . 
x 


(A-Z) 


If  oa  =  0, 

y 


b  =  0 

y 


0  5  b  <:  1 
x 


(A-  3) 


b  =  0 

y 


-pa  I 

1  ^  b  ^ 

X  I  —  I 
x  y 


(A  -4) 


If  na  >  0  and  pa  >  0, 
x  y 


pa  I  b 

b  =  y  x  x 

y  "pa  I  +  (I  —I  )b 
x  y  x  y  x 


0  <  b  <1 
x 


(A-5) 


pa  I  b 
y  x  x 

V  "  pa  I  +  (I  -I  )b 
x  y  x  y  x 


a  I  \l/3 
x  y  \ 


-  1  <  b  s  -b  4 

x  y.  a  I 

\  y  x. 


-oa  I 

and  b  ^  T  X/ 

X  I  —  1 

x  y 


(A-6) 


These  expressions  denote  surfaces  in  b  ,  b  ,  b  space.  On  the  b  ,  b  plane 

x  y  z  x  y 

they  denote  curves  which  are  shown  in  Figs.  A-l,  A-Z,  A-3,  and  A-4.  To 
derive  these  expressions,  let 


T(b  ,  b  )  =  T./h1  +  T,/hZ 

x  y  1  Z 


where 


T.  =  (b  -pa  )2/l  +  (b  -pa  )Z/l 

i  xxx  y  y  y 


13 


and 


T2  =  -  0*//lz+02/C. 


Define  r  and  0  such  that  b  =  r  cos  0  and  b  =  r  sin  0,  and  consider 

x  y 

circles  of  constant  r  on  the  b^,  b^  plane.  A  necessary  condition  for  a  point 

at  radius  r  to  be  a  minimum  of  T(b  ,  b  )  is  that  on  a  circle  of  constant  r 

x  y 

this  point  be  a  minimum  of  T(0),  a  function  of  a  single  variable  0.  Since 
is  independent  of  0,  a  point  is  a  minimum  of  T(0)  if  and  only  if  it  is  a  mini¬ 
mum  of  Tj(0).  The  conditions 


and 


1  dTi 

2  W 


1  d'Tl 


b  (b  -pa  ) 
x  y  y 


b  (b  -pa  ) 
y  x  x 


=  0 


x 


■  b  (b  -pa  )  +  (b  ) 

y  y  y  x 


2  de2 


b  (b  -pa  )  —  (b  ) 
XX  x  y 

I 

x 


*  0 


which  must  be  satisfied  by  all  minima  for  r  >  0  yield  the  desired  expressions. 
The  point  r  =  0  must  be  checked  separately. 


2.  All  minima  of  T(b  ,b  ,b  )  satisfy  the  following  expressions.  If 

x  y  z 


Dax  =  °' 


b  =  0 
x 


pa  I 

o  <;  b  ^  r^r 

z  I  —  I 
X  z 


(A-  7) 


Daz  Xx 

n—=TT 

X  z 


-  1  <  b  5  1 
x 


(A  -  8) 


If  pa  =0, 
z 


=  0 


0  <  b  si 
x 


(A-9) 


=  0 


-pa  I 

,  ,  x  z 

1  s  b  s  -j - =-- 

x  I  —  I 
x  z 


(A-  1  0) 


If  pa  >0  and  pa  >0, 
x  z 
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(A  -  1  1 ) 


b 

z 


pa„ 


I  b 
x  x 


pa  I  +  (I  —I  )b 

X  Z  X  Z  X 


0  <  b  <  1 

X 


b 

z 


pa  I  b 

Z  X  X 

pa  I  +  (I  —  I  )b 

X  Z  X  Z  X 


—  1  <  b  < 

X 


1/3 


and  b  < 
x 


-pa  I 
x  z 

I  -I 

X  z 


(A  -  1  Z ) 


These  expressions  are  identical  to  those  of  Section  1  except  that  y  is  replaced 
by  z.  They  are  derived  in  the  same  way  as  the  expression  of  Section  1. 


3.  Combining  the  results  of  Sections  1  and  Z  yields  curves  in  three- 

dimensional  space  on  which  all  minima  of  T  must  lie.  These  curves  are 

inte  r  sections  of  the  surfaces  of  Sections  1  and  Z  and  are  shown  in  Figs.  A-5 

through  A-10  for  six  cases  that  occur.  A  seventh  case,  a  =  0,  a  >0,  and 

x  y 

a^  >  0  which  is  shown  in  Fig.  A- 11  requires  further  study.  The  results  of 
Sections  1  and  Z  imply  that  any  minimum  of  T  must  be  either  on  the  straight 
line 

^  oa  I 

0:b25l,  b  = 

x  y  I  —  I 

3  x  y 

or  the  rectangle 

pa  I 

b  =  0,  0  <;  b  g  T  y  T- 

x  y  I  —I 

3  x  y 

Differentiating  the  expression  for  T  on  this  rectangle  yields  further  results 
that  all  minimum  points  on  the  rectangle  satisfy 


b  = 


pa  I 
z  x 


z  I  —  I 
x  z 


pa  I 

0  <  b  <  -  — 

z  I  —  I 
X  z 


b 

y 


pa  I  b 
y  z  z 

I  +  (I  -I  )b 


Pa. 


y 


y 


0  <  b 

z 


pa  I 
z  x 


(A- 13) 


This  curve  is  also  shown  in  Fig.  A- 11. 

Z  Z  Z 

The  constraint  b  +  b  +  b  =1  implies  that  only  values  of  T  defined 
x  y  z 

for  points  on  the  sphere  of  unit  radius  centered  at  the  origin  are  of  interest. 

In  particular,  all  minimum  points  of  T  are  on  this  unit  sphere  so  that  for  each 


1  5 


case  of  Figs.  A-4  to  A-ll  a  minimum  of  T  can  only  occur  where  a  curve 

intersects  the  unit  sphere.  Therefore,  there  can  be  at  most  two  or  three 

minima  of  T  depending  on  the  case.  There  is  always  one  minimum  in  the 

region  ^  0,  b^  ^  0,  b^  ^  0  and  there  can  be  zero,  one,  or  two  minima  in 

the  region  b^  <0,  b^  >  0,  b^  >  0.  Equations  A-l  through  A- 13  may  be  used 

to  determine  the  coordinates  of  all  minima  of  T.  The  next  section  investigates 

sufficient  conditions  for  there  to  be  only  one  minimum  or  equivalently  there  to 

be  no  minima  in  the  region  b  <  0,  b  >  0,  b  >0. 

*  x  y  z 

4.  From  Eqs.  A-l  through  A-6  or  Figs.  A-l  through  A-4,  it  is  evi¬ 
dent  that  sufficient  conditions  that  no  minima  occur  in  the  region  b  <0, 

b  ^  0,  b  2?  0  are  as  follows, 

y  z 

If  pa^  =  0, 


-oa  I 

x  y 

I  -I 

x  y 


-  1. 


If  pa^  =  0, 


pa 

y 

m 

X 


I 

X 


I 

y 


< 


l. 


If  pa  >  0,  pa  >0, 
x  y 


(d  .  r  >  i 

mm 


where  d  .  is  the  minimum  value  of 
min 


d  = 


pa  I  b  \ 

2 

b  Z  4 

X 

_ 

y  x  x 

oa  I  +  (1  —  I  )b 

x  y  x  y  x 

1/2 


the  distance  from  the  origin  to  a  point  on  the  surface  on  which  all  minima  must 

2 

occur  in  the  region  b  <  0,  b  ^  0,  b  ^  0.  Differentiation  of  d  with  respect 

x  y  z 

to  b  yields 
x 

,  ,  (a  I  )2/3  +  (a  I  )2/3'  3 

2  2  x  y _  y  x 


(d  .  )  =  o 

min 


(I  -I  y 

x  y 
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In  all  three  cases,  the  sufficient  condition  for  no  minima  in  the  region  b  <  0, 
b  5  0,  b  5  0,  which  implies  that  there  is  only  one  minimum  of  T,  is 

y  z 


Z 

P  > 


<1  -I  ) 

je y 


(a  I  )2/3  +  (a  I  )2/3]3 
x  y  y  x 


An  analogous  sufficient  condition  for  there  to  be  only  one  minimum  of  T  can 
be  derived  from  Eqs.  A-7  through  A-1Z.  It  is 


o2  > 


(I  -I  ) 

x  z 


[la  I  )2/3  +  (a  I  )2/3"|3 
[_  X  Z  Z  X  J 
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P°t~  o 


3-66-7184 


P  °y  =  0 
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P°x  >0,  p oy  >  0  ,  />> 


( a x I y )2/3  (oyIx)1/3 


3-66-7185 


P(sx> 0,  /5ay>0,  />< 


(axIy)2/3  (ayIx)'/3 


3  66-71  86 
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ax  =  1 ,  ay  =  a2  =  0 


3-66-7187 


b2 


3-66-7188 


ay  =  1 .  az  =  0 
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3-66-7189 


az  =  1 ,  ax  =  ay  =  0 


3-66-7190 


Oy^O,  ax  >0,  a2  >0 


b* 
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oz  =  0,  ax >0,  az>0 


3-66-7191 


P°  ylx/dx- 
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ox>0,  oy>0,  az>0 


3-66-7192 
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3-66-7193 


Ox  =  0,  0y >0, oz  >0 


b. 


Fig.  A  - 1 1 
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